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迭代 型 去 模糊 算法 的 自 适应 迭代 终止 条 件 
江 顺 亮 ， 姜 尹 楠 ， 曾 小 霞 ， 唐 祥 玲 ， 徐 少 平 ; 


(南昌 大 学 信息 工程 学 院 , 南昌 330031) 


摘 X: 由 于 缺乏 有 效 的 迭代 终止 条 件 (iterative stopping criterion，ISC)， 目 前 大 多 数 去 模糊 算法 简单 采用 固定 的 迭代 
次 数 实现 ， 存 在 着 执行 效率 低 、 去 模糊 效果 不 够 理想 等 问题 。 为 此 ， 提 出 一 种 基于 残 差 图 像 ( 选 代 过 程 中 Home 中 间 
估计 图 像 和 模糊 核 卷 积 后 与 模糊 图 像 之 间 的 差 值 ) 的 去 模糊 效果 度量 (deblurring measure, DM), 并 在 该 DM 度量 的 基础 
上 设计 了 一 种 自 适应 的 迭代 终止 条 件 (adaptive ISC，AISC) 。 将 所 提出 的 AISC 迁 代 终止 条 件 应 用 于 经 典 的 
NCSR(nonlocally centralized sparse representation) 迭代 型 去 模糊 算法 中 。 在 均匀 模糊 、 高 斯 模糊 和 运动 模糊 三 种 典型 模 
糊 失 真 类 型 下 大 量 的 实验 数据 表明 ， 与 采用 国定 和 迭代 次 数 的 原 NCSR 算法 相 比 ,采用 自 适 应 迭代 条 件 后 NCSR 算法 执 
行 效率 得 到 显著 提升 ， 且 所 复原 图 像 在 PSNR、SSIM 和 FSIM 图 像 指标 值 上 与 原 算 法 差别 不 大 。 
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Adaptive iterative stopping criterion for deblurring algorithm 


Jiang Shunliang, Jiang Yinnan, Zeng Xiaoxia, Tang Yilingm, Xu Shaoping’ 
(School of Information Engineering, Nanchang University, Nanchang 330031, China) 


Abstract: Currently, due to lack of effective iterative stopping criterion (ISC) , most of the ID algorithms are usually 
implemented with fixed iteration steps, leading to the low efficiency and poor deblurring results. Therefore, this paper proposed 
a deblurring measure (DM) based on the residual image (the difference between the intermediate estimated image convoluted 
with blur kernel and blur image during each iterative step) . Based on the proposed DM metric, it designed and applied an 
adaptive iterative stopping criterion (AISC) to the classic nonlocally centralized sparse representation (NCSR) . Extensive 
experiments on uniform blur, Gaussian blur and motion blur show that, compared with the original NCSR algorithm using fixed 
iteration steps, the NCSR algorithm adopting AISC has an obvious advantages in terms of efficiency and the 
reconstructed images have a similar image quality in terms of peak signal-to-noise ratio (PSNR) , structural similarity (SSIM) 
and feature similarity (FSIM) . 
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0 引言 欠 约 束 条 件 问题 ， 得 到 了 众多 学 者 的 关注 并 已 经 提出 了 很 多 图 
像 去 模糊 (image deblurring) 算 法 上 09。 其 中 ， 基 于 正则 化 的 求解 
图 像 模糊 失真 过 程 一 般 可 以 形式 化 的 表示 为 方法 是 一 类 比较 常见 和 有 效 的 解决 方案 。 在 基于 正则 化 的 去 模 
g-f*h Q) WEH, SCEDRUHUSUG EE. f 的 先 验 知识 构造 合适 的 正则 项 ， 
Mp: g 和 分 别 是 模糊 图 像 ( 即 观测 图 像 ) 和 原始 图 像 ， 有 是 。” 并 对 求解 的 复原 图 像 和 退化 函数 进行 正则 化 约束 。 其 数学 模型 
点 扩散 函数 (point spread function,PSF)， 通 常 称 为 退化 函数 ，* DAWR O2) 所 示 。 
为 卷 积 运算 符 。 图 像 去 模糊 是 一 个 典型 的 求 逆 问题 ， 主 要 基于 1 
i a S ue : C,(f.h) - -lg - f * n]; + A99) Q 
观测 到 的 模糊 图 像 通过 求解 关于 某 个 目标 函数 的 最 优 值 实现 对 2 
原始 无 失真 图 像 的 最 佳 估计 。 由 于 图 像 去 模糊 问题 属于 典型 的 ”其 中 :第 一 项 为 数据 保 真 项 , eo) 是 正则 项 表示 图 像 先 验 知识 ; 
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E 则 化 参数 。 求 解 式 (2 站 
variant, TV) [5 一 类 方法 因 具 备 良 好 的 
能 力 ,得 到 了 广泛 的 应 月 


日 。 但 由 于 忽略 了 


FP 所 定义 的 


优化 问题 ， 全 变 分 (total 
持 和 抑制 噪声 的 


细节 与 结构 特性 ， 


means) 降 


算法 ， 解 决 了 基 
中 ， 利 用 图 像 的 非 局 部 
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TV 去 模糊 
REATHA R, V 
引入 到 正则 化 方法 中 ， 提 出 了 基于 


容易 造成 图 像 平坦 区 域 产 生 阶梯 效应 (staircasing effect), 所 以 许 
FE 来 , 受 NLM (non-local 
局 部 自 相 似 性 


算法 被 提出 来 。 近 4 
F 多 学 者 把 
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似 性 
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(normalized graph laplacian) 


E 则 化 约束 的 去 模糊 
算法 中 的 不 足 &5?1。 在 这 些 


的 图 像 局 部 空 


或 ， 充 分 利用 了 图 像 


很 大 的 改善 。 


$ A PIAJ H JH 


以 的 特点 。 该 算法 通过 


gon Ne fe E 


作为 正则 
内 外 两 层 和 迭代 过 


其 于 先前 估计 来 更 新 村 


计算 相似 性 


项 ， 充 分 发 挥 了 空间 自 相 
F 程 实现 ,在 外 部 迭代 中 ， 
E(similarity weights), 内 


部 迭代 使 用 


图 像 的 非 局 部 
原 算法 0o-331。 例 如 ，Zhang 
的 稀疏 表示 在 算法 执行 效率 和 稀 玻 


Jeu pe RETE E rds] NR 


H AHATE- Bib PE FH 


提出 许多 新 的 图 像 
A 传统 基于 图 像 块 


等 人 0 为 了 解 y 


标 函 数 。 最 近 ， 学 者 们 


于 以 


模糊 化 效果 。 
b) 去 模糊 后 的 图 像 质 


较 少 迭代 步骤 完成 后 ， 和 迭代 中 间 过 程 产 4 
量 值 就 达到 了 最 佳 。 然 后 随 着 迭代 次 数 的 增 力 
而 下 降 。 
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值 仅 是 缓慢 上 升 ， 直 到 720 次 迭代 完成 所 得 到 的 图 像 质量 值 
提升 量 很 小 (视觉 上 不 易 感觉 到 差别 ) 。 换 句 话 来 说 ， 采 
和 迭代 次 数 的 方法 浪费 了 较 多 的 执行 时 间 ， 却 没有 大 幅度 提高 去 


未 必 是 最 佳 的 。 在 某 


些 时 候 ， 经 过 


过 程 中 利用 


ES UL. fu 


H 
ZR 


般 比 较 高 。 


忆 此 ， 上 述 两 种 情况 表明 需要 和 
NCSR 算法 这 类 人 迭 代 型 去 模 糙 
以 少 的 执行 代价 获得 最 佳 的 去 模 狗 
解决 NCSR zz bl 


前 评价 效果 


算法 上 述 问 
无 参考 的 图 像 质量 讨 
quality assessment, NR-IQA) 对 每 一 步 迭 代 产 生 的 中 间 估 计 图 像 


究 一 种 迭代 终止 条 件 应 用 
算法 中 ， 使 得 这 些 算法 能 够 
化 效果 。 


题 一 个 简 


甘于 组 稀疏 表示 的 (group sparse reprentation, GSP) 图 像 复 


GSP(group sparse regularization) 
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可 像 去 模糊 上 。 该 算法 主要 利用 
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与 非 局 部 自 相 似 结合 ， 


用 Bregman 4) £85 [V (split bregman iteration, SBD) 方 法 求解 复 
i Ee ER P e o ips 类 似 上 


H, Wang 等 人 0 提出 基于 迭代 增 量 的 


式 更 新 对 原始 图 像 的 估计 ; Liu 等 人 
忧 合 稀疏 表示 (joint sparse reprentation, JSR) 算 法 。JSR 通过 


正则 化 去 模糊 算 
非 负 约束 和 稀 
CHE 3D 转换 域 上 提出 了 


大 例 , 算法 采 月 


法 , 运用 在 遥感 
iL ZW EET 


实现 获得 最 佳 的 去 模糊 图 像 的 质 
的 执行 效率 同时 急剧 下 降 。 
AGE Ré rp Br nn gap eei E 
昌 在 迭代 过 程 中 中 间 估 计 医 


为 此 ， 本 文 不 采用 


BA 


NR-IQA 算法 


E 的 中 间 估 计 图 像 的 质 
1， 图 像 质量 值 反 


方法 就 是 在 迭代 


F 价 算法 (no-reference image 


进行 评价 ， 根 据 评价 值 的 高 低 决定 迭代 过 程 是 否 可 以 终止 。 但 
于 NR-IQA 算法 并 非 专门 为 去 模糊 算法 而 设计 ， 其 计算 复杂 
比较 好 的 NR-IQA 算法 本 身 的 执 
行 时 间 比 较 长 , 大 多 达到 秒 的 量 级 , 有 些 甚 至 高 达 几 十 秒 05-17。 
在 NCSR 算法 和 儿 代 求解 过 程 中 反复 频繁 调用 


将 会 导致 去 模糊 算法 本 身 


cr 


的 差 值 图 像 ( 本 文 将 差 人 
的 迭代 终止 条 们 
也 ， 基 于 残 差 
于 去 模糊 算法 上 


图 像 中 像素 点 之 间 亮 度 


以 去 模糊 算法 迭 
F 价 对 象 的 常规 做 法 ， 而 是 
像 和 模糊 核 卷 积 后 与 模糊 图 像 之 
直 图 像 称 为 残 差 图 像 ) 设 计 一 种 自 适 应 


F(adaptive iterative stopping criterion, AISC) 。 


aj ACE Em. (EAE 


RAI RA EE ERKAK 
执行 效率 不 高 和 去 模糊 效果 未 必 最 佳 的 问题 。 
目 了 内 外 两 层 循环 , 外 循 
层 循环 为 120 次 ， 共 有 720 次 迭代 。 
NCSR 算法 采用 固定 的 迭代 次 数 的 实 


下 就 已 达到 了 比较 高 的 


岗 方 法 存在 着 


以 NCSR 算法 [1 
迭代 次 数 为 6 次 ， 
通过 实验 可 以 发 现 ， 
岗 方法 会 导致 以 下 两 种 情 


上 对 候 ， 该 算法 在 较 少 的 迭代 次 数 
图 像 质量 值 。 在 后 继 的 迭代 中 ， 图 像 质 


对 每 个 图 


线 构造 字典 ， 然 后 将 模糊 图 
Bolt fT LARES. A 


值 差 异 大 小 定义 了 专用 
评价 去 模糊 效果 的 度量 (deblurring measure, 


非 局 部 的 规则 项 ， 能 有 效 刻 画图 像 的 非 局 部 统计 特征 ， 相 DM), 然后 基于 该 DM 度量 设计 了 AISC HEMARA ERF. 
对 现 有 的 稀 玻 表示 方法 更 为 准确 ; Liu 算法 在 贝 叶 斯 框架 下 使 为 验证 AISC 的 实际 效果 ， 将 所 提出 的 AISC RRA IEAR TIEN 
I 迭代 方法 来 求解 关于 原始 未 失真 图 像 的 最 优 解 ,Dong 等 有 于 经 典 的 NCSR 算法 中 。 在 均匀 模糊 、 高 斯 模糊 和 运动 模 站 
是 出 的 NCSR(nonlocally centralized sparse representation) $ HE 型 模糊 失真 类 型 下 大 量 师表 明 ， 与 采用 固定 迭代 次 
法 把 关于 图 像 稀疏 先 验 知识 与 非 局 部 相 E 一 起 ， 首 先 利 数 的 原 算法 相 比 ， 采 用 AISC ARZE NCSR 算法 所 复原 图 
辐 像 非 局 部 重复 的 特性 估计 噪声 图 块 的 稀 玻 系数 ， 然 后 采用 像 在 PSNR、SSIM08 和 FSIM09 图 像 指标 值 上 比 原 算 法 相差 不 
性 缩 图 块 稀 朴 系 数 的 方法 实现 图 像 的 复原 ， 相 比 与 大 多 数 大 ， 而 执行 效率 则 显著 得 到 提升 。 
算法 都 有 更 好 的 复原 效果 。 
目前 ， 就 基于 正则 化 技术 实现 的 各 类 去 模糊 算法 而 言 ， 通 相关 工作 
常 通过 迭代 方法 来 求解 目标 函数 最 优 图 像 的 最 佳 估 NCSR 算法 简介 
计 [ 鸣 。 然 而 由 于 缺乏 判断 迭代 终止 的 条 件 ， 这 些 算法 常常 只 能 在 经 典 的 NCSRD3] 算 法 中 ， 先 基于 给 定 的 模糊 图 像 本 身 在 


像 分 解 成 图 块 并 基于 所 学 习 的 字典 


收缩 (shrinkage) 方 法 处 理 每 个 


定 字典 情况 下 计算 出 
f 2 Rug 表示 在 像素 点 ! Able 
ÆA g 中 像素 点 i 处 
peR”” n<M> MA 


RARR SUELE, P 


实现 图 像 去 模糊 。 其 中 ， 


[| 


P PALM EL ERU. XF 


了 按照 上 述 方法 的 逆 过 程 重 构图 像 从 而 
稀 琉 编码 是 核心 技术 ， 其 任务 是 在 给 


RADY Jax Vi 的 
ERES y 的 矩阵 操作 。 给 定 字 典 
图 像 中 提取 的 某 个 图 块 可 稀疏 表示 为 


图 像 。 RYG 


图 块 , HE R 
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录用 稿 


g = GQ,;， 图 像 的 去 模糊 问题 通过 求解 


min, > ( 
其 中 : p 是 从 图 像 g 中 提取 的 图 块 总 数 ， w, 相对 于 字典 jg 而 
言 是 指 图 块 g, WARRI: A 表示 正则 化 参数 。 图 像 g 可 以 


如 下 的 优化 问题 实现 


G) 


lg; - $a, à 


通过 稀疏 系数 集合 {a,,) 来 表示 。 从 集合 [a } 对 图 像 8 进行 重 
建 是 一 个 超 定 系 统 (over-determined system) 问 题 , 可 采用 最 小 二 
乘 方法 求解 : 
f ~ RERS (RGa,,) a 
为 便于 表示 ， 式 (5) 可 改写 为 
g~po0, =Q (Rta) (5) 
IP: a, 表示 所 有 w ,的 级 联 。 这 样式 (4) 可 以 写 为 
a, -argmin, (|g; -pe° al, + 4|a]) ©) 
在 NCSR 算法 中 , 通过 解决 下 面 的 最 小 化 问题 实现 从 g fih 
i f: 


Q; 三 arg min, (|; -go C 


CAM lo; -Bl,) (7) 
其 中 : 4 是 平衡 保 真 度 与 稀疏 项 的 参数 ， p RERO TARE 
数 w 的 估计 。NCSR 算法 应 用 于 图 像 去 模糊 ， 主 要 依靠 迭代 过 
程 求解 式 (8)。 在 每 一 步 迭 代步 又 中 都 会 产生 关于 原 图 像 x 的 估 
计 值 


jet = ÊO eh (g - hf?) (8) 


其 中 :5 是 预定 义 的 常量 ， 8 是 模糊 图 像 ， 是 模糊 核 ， po 


用 于 获 


是 第 了 次 适 代 所 估计 的 图 像 ， PORRER, f 


得 j+1 步 的 稀 玻 编码。 去 模糊 主要 依靠 收缩 方法 计算 新 的 a1 
实现 ， 定 义 如 下 : 


aj" 2 S,(a] - B) + B (9) 

式 (10) 中 的 s, C) ZEE PEERS s c 是 软 阔 值 参数 。 在 每 

次 内 层 迭 代 中 ，w 六 一 旦 获得， 去 模糊 的 中 间 图 像 就 可 以 重 
为 


a 


f0 = poat” (10) 
需要 说 明 的 是 ， 在 NCSR 算法 中 ， 固 定 地 采用 pon 作为 

原 图 像 的 最 终 估 计 值 ( 即 去 模糊 后 的 图 像 )。 
12 存在 问题 
在 和 迭代 次 数 的 选择 上 ，NCSR 算法 是 以 设置 固定 值 的 方式 
来 处 理 的 。 然 而 各 种 图 像 的 内 容 和 模糊 失真 类 型 有 着 各 自 的 特 
点 ， 在 算法 中 国定 设置 的 720 次 迭代 不 能 保证 对 不 同 内 容 的 图 
像 和 各 类 失真 都 能 够 得 到 很 好 的 去 模糊 效果 。 为 了 分 析 NCSR 
算法 采用 固定 迭代 次 数 实现 所 产生 的 问题 ， 采 用 常用 测试 图 像 
Barbara( 大 小 为 256*256) 为 例 , 对 算法 在 均匀 模糊 、 高 斯 模糊 、 
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运动 模糊 三 种 典型 模糊 失真 下 在 720 次 迭代 过 程 中 图 像 的 图 像 
RAE PSNR, SSIM, FSM 的 变化 趋势 。 实 验 的 参数 设置 中 ， 

保证 模糊 类 型 和 失真 程度 大 小 的 全 面 性 。 在 均匀 模糊 下 ， 
模糊 核 (Kernel size, 简 记 K) 大 小 设置 为 [5 5]、 [15 15]、 [2525]: 
高 斯 模糊 参数 设置 为 :模糊 核 大 小 K 设置 为 [5 5]. [9 9]. [15 
15], 标准 偏差 Sigma( 简 记 S) 为 1 和 3; 运动 模糊 参数 为 : 距离 
Length( 简 记 也 为 10、20、40， 角 度 Angel( 简 记 A) 为 0。 表 1 数 
据 记 录 了 三 种 模糊 类 型 下 ，NCSR 算法 在 某 些 关键 迭 代步 骤 上 
的 图 像 质量 值 。 从 表 1 数据 可 以 看 出 ， 图 像 质量 值 有 两 种 变化 
趋势 : DW A， 图 像 质量 值 在 前 面 较 少 的 迭代 次 数 时 便 达 到 
峰值 , 在 随后 的 迭代 次 数 中 , 图 像 质量 值 一 直下 降 ; b) 情 况 B， 

图 像 质量 在 迭代 前 期 提升 很 快 ， 在 后 期 的 迭代 次 数 非常 缓慢 的 
上 升 。 对 于 情况 A 来 说 ， 算 法 发 生 了 过 和 迭代 的 现象 ， 过 多 的 迁 
代 不 仅 给 算法 造成 了 很 大 的 时 间 开 销 ， 而 且 随 着 和 迭代 次 数 的 增 
加 复原 的 图 像 质量 越 来 越 低 ， 就 情况 B 而 言 ， 算 法 在 前 面 的 迭 
代 中 图 像 质量 值 变化 很 大 ， 有 非常 明显 的 提升 ， 而 后 期 次 迭代 
中 图 像 质量 值 变 化 缓慢 , 图 像 质量 值 变 化 不 大 (从 人 眼 视 觉 观察 
角度 )。 总 之 ,无论 是 对 于 上 述 哪 种 情况 而 言 ， 需 要 自 适 应 的 设 
迭代 终止 条 件 下 从 而 保障 NCSR 算法 在 获得 最 佳 的 去 模糊 图 
像 质量 的 同时 具有 令 人 满意 的 执行 效率 。 

为 了 更 为 直观 地 反映 复原 图 像 随 着 迭代 变化 的 趋势 ， 图 1 
分 别 列 出 了 两 种 情况 下 Barbara 图 像 的 去 模糊 效果 。 其 中 图 1 
(a) (b) 反 映 的 是 上 述 变 化 趋势 中 的 A 情况 。 从 图 (a) ( Gs 
动 模糊 A=0, L=10) 可 以 发 现 , 720 次 迭代 去 模糊 后 的 图 像 ( 左 图 ) 
明显 较 最 佳 质 量 值 图 像 ( 右 图 ) 要 模糊 ; 图 1(c)(d) 反映 的 是 上 述 
变化 趋势 中 的 B 情况 ， 在 图 1(c)(d) (高 斯 模糊 K-[15 15], S=3) 
中 , 左 图 为 720 次 迭代 后 获得 的 图 像 , 右 图 为 整个 去 模糊 过 程 中 
最 佳 的 图 像 (在 迭代 200 次 获得 )。 可 以 发 现 ， 两 张 图 像 的 差别 


E 
FE 
E 
FH 


非常 小 ， 任 人 眼 无 法 辨别 出 差别 。 由 此 说 明 ，NCSR 算法 需要 
添加 一 个 迭代 终止 条 件 ， 在 合理 的 友 代 次 数 下 实现 完成 效率 和 


复原 图 像 质量 之 间 的 平衡 。 
2 ” 自 适 应 迭代 终止 条 件 


2. RRI 
由 上 文 可 知 ,本 文 所 指 的 残 差 图 像 是 指 在 途 代 过 程 中 的 中 
间 估 计 图 像 和 模糊 核 卷 积 后 与 模糊 图 像 之 间 的 差 值 图 像 。 
r2g-f*h 1) 
理论 上 讲 ， 去 模糊 算法 在 刚 开 始 进行 迭代 的 过 程 中 ， 由 于 
对 于 原 图 像 的 估计 不 是 很 准确 ， 这 个 时 候 残 差 图 像 中 各 个 像 
点 之 间 在 差异 比较 大 (尤其 是 较 难 复原 的 图 像 边缘 细节 部 分 )。 
a 至 渐 增加 , 估计 图 像 # 越 来 越 准确 , 则 Pap 越 
来 越 与 失真 图 像 g 近似 ， 这 样 导致 残 差 图 像 p 图 像 中 各 个 像素 
点 亮度 值 变 小 且 趋 于 在 0 值 附近 变动 ,像素 点 亮度 值 之 间 的 差 
异 也 趋 于 变 小 。 以 一 张 256*256 大 小 的 Barbara 图 像 为 例 ， 分 
别提 取 去 模糊 迭代 过 程 中 第 1、10、20、50 和 720 共 5 个 迭代 
位 置 的 残 差 图 像 ，5 个 迭代 位 置 的 残 差 图 像 如 图 2 所 示 ( 为 便于 
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IH 


视觉 观察 ， 亮 度 值 已 经 规范 化 到 [0 255] 间 。 限 于 篇 幅 ， 这 里 只 差 图 像 像素 点 亮度 值 之 间 的 差异 越 来 越 小 在 选 代 次 数 为 50 时 ， 
给 出 运动 模糊 类 型 下 (A=0,L=10) 的 残 差 图 像 图 ， 均 匀 模 糊 与 高 。” 从 视觉 上 已 经 无 法 观察 出 残 差 图 像 中 有 关于 图 像 内 容 的 信息 ， 
斯 模糊 情况 类 似 )。 说 明 图 像 去 模糊 效果 接近 最 优 值 。 图 3 是 上 述 5 个 迭代 步 又 
图 2 可 以 看 出 ， 在 和 迭代 刚 开始 的 第 一 步 时 ， 残 差 图 像 中 残 差 图 像 像素 点 亮度 值 对 应 的 直方 图 。 由 图 3 可 知 ， 残 差 图 像 
明显 看 到 能 反映 图 像 内 容 的 关键 边缘 信息 。 到 第 20 步 的 时 候 ， ”中 亮度 值 的 变化 范围 随 着 迭代 次 数 增加 ， 在 逐渐 收 窗 ， 说 明 中 
尚 能 看 到 不 明显 的 结构 信息 。 之 后 ， 随 着 迭代 次 数 的 增加 ， 残 。 间 估 计 图 像 在 逐渐 接近 原 图 像 。 
表 1 Barbara 图 像 去 模糊 过 程 中 图 像 质量 值 变化 趋势 
模糊 类 型 。 参数 。 ”评价 标准 ”第 1 步 第 10 步 第 20 步 第 50 步 第 720 步 变化 类 型 


pmi 


p 


PSNR 22.65 29.42 29.60 29.14 28.70 A 
HB) K=[5 5] SSIM 0.6127 — 0.8433 0.8303 0.8035 0.7856 A 
FSIM 0.7472 | 0.9057 0.9243 0.9152 0.9075 A 
PSNR 20.84 23.11 23.52 23.91 24.64 B 
均匀 模糊  K-[1515] SSIM 0.4554 © 0.6065 0.6296 0.6510 0.6878 B 
FSIM 0.6173 0.7603 0.7868 0.8082 0.8349 B 
PSNR 19.30 21.72 22.14 22.64 23.67 A 
均匀 模糊 ”K=[25 25] SSIM 0.3871 0.5156 0.5387 0.5644 0.6211 B 
FSIM 0.5568 0.6696 0.6989 0.7270 0.7612 B 
PSNR 22.99 28.70 29.01 28.18 25.98 A 

K=[5 5] 
高 斯 模糊 SSIM 0.6369 0.8292 0.8081 0.7612 0.6763 A 
2i FSIM 0.8005 | 0.9177 09141 0.8969 0.8622 A 
PSNR 22.84 24.32 24.79 25.38 25.81 B 

K=[9 9] 
高 斯 模糊 SSIM 0.5909 0.6877 0.7029 0.7072 0.6796 A 
FSIM 0.7216 0.8165 0.8355 0.8527 0.8522 A 
PSNR 22.30 23.46 23.55 23.66 23.89 B 

K=[15 15] 
高 斯 模糊 SSIM 0.5571 0.6416 0.6483 0.6546 0.6639 B 
d FSIM 0.6886 0.7829 0.7908 0.7981 0.8087 B 
PSNR 23.02 26.61 27.12 26.66 25.38 A 
运动 模糊 SSIM 0.5943 0.7640 0.7486 0.7053 0.6560 A 
s FSIM 0.7655 0.8771 0.8780 . 0.8645 0.8448 A 
PSNR 20.99 23.77 24.18 24.41 24.75 B 
运动 模糊 = SSIM 0.4751 0.6485 0.6496 0.6293 0.6263 A 
ic FSIM 0.6806 — 0.8130 — 0.8231 0.8211 0.8201 A 
PSNR 18.92 22.48 23.17 23.55 24.03 B 
运动 模糊 i SSIM 0.3944 — 0.5803 0.6079 0.5980 0.5917 A 
Mi FSIM 0.6105 — 0.7626 0.7920 0.8038 0.8069 B 


(d) 高 斯 模糊 (K=[15 15], S-3) F 


(a) 运动 模糊 (A=0, L-10)F ——— (b) 运动 模糊 (A=0, L-10) F 


(c) 高 斯 模糊 (K=[15 15], S23) 
第 720 步 对 应 的 图 像 最 佳 质量 值 对 应 的 图 像 下 第 720 步 对 应 的 图 像 ”第 200 步 获得 的 最 佳 图 像 质量 的 图 像 
图 1 迭代 过 程 中 不 同 迭 代步 又 下 所 获得 的 图 像 质量 在 视觉 上 的 对 比 
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(a P HUS 
1 时 残 差 图 像 


图 


(b ACERO (c) P HUS 
10 时 残 差 图 像 20 时 残 差 图 像 50 时 残 差 图 像 
图 2 去 模糊 时 关键 迭代 步 数 下 所 得 到 的 残 差 图 像 


120r 


(MAREKO 


(REA 
720 时 残 差 图 像 


100 上 1 100 
80r 1 80 
60r 1 60 
40} 1 40 
20! 1 20 
0 o! IB.. 0 m" bi.. 
-01 -008 -006 -004 -002 0 002 004 006 008 -003 这 ül 00? 00 
人 图 像 直方 图 m 方 图 


120r 


d 2 M li 4005 0 0005 bi 0015 002 0.025 


(dA REC 50 时 残 差 
3 ”去 模糊 时 关键 迭代 步 数 下 所 得 到 


CHEIRA 10 时 下 


图 


0 
-002 -0015 -001 -0005 0 0.005 001 0.015 002 0.025 


像 直方 图 图 


(e) 迭 代步 数 为 720 FERA: 


图 


f£ El 


图 


方 


的 残 差 图 像 的 直方 图 


2.2 和 迭代 终止 度量 由 式 (13) 的 定义 可 知 ，Coe (m,n) 的 值 越 小 , 表明 像素 点 2 
基于 图 2 和 3 对 残 差 图 像 中 各 个 像素 点 亮度 值 的 变化 趋势 ” 间 的 差异 越 小 ,为 了 计算 邻 域 范围 内 像素 点 亮度 值 的 差异 程度 ， 
分 析 可 以 知道 ， 随 着 迭代 次 数 的 增加 ， 去 模糊 后 的 估计 图 像 越 。” 可 以 定义 QL+1)*(2L+1) 邻 域 窗口 范围 下 中 心 像素 点 与 邻 域 内 
来 越 接近 原 图 像 ， 这 样 导 致 残 差 图 像 中 各 个 像素 点 的 亮度 值 变 ”各 个 像素 点 的 差异 总 和 DM ,作为 衡量 复原 图 像 整体 去 模糊 效 
小 ， 之 间 的 差异 程度 也 越 来 越 小 。 基 于 此 ， 本 文采 用 下 面 的 方 果 的 度量 。 
法 定义 了 一 种 衡量 去 模糊 效果 好 坏 的 DM 度量 。 具 体 地 ,根据 " 2 
式 (11) 计 算出 残 差 图 像 , 然后 对 残 差 图 像 进行 归 一 化 , 归 一 化 后 0 (14) 
残 差 图 像 的 均值 为 零 ， 方 差 为 1。 MR 
"NE (12) 其 中 : 工 是 邻 域 范围 的 半径 。 表 2 是 大 小 为 256*256 的 Barbara 
Varr) 图 像 分 别 在 均匀 模糊 (K=[5 5])、 高 斯 模糊 (K=[15 15], S=3) 和 运 
Hop: /代表 归 一 化 后 的 残 差 图 像 ， 地 和 var(r) 分 别 是 残 差 图 动 模糊 (A=0, L=10) 中 的 度量 值 DM ,的 变化 情况 。 从 表 2 中 数 
像 的 均值 和 残 差 图 像 的 方差 。 基 于 此 ， 假 设 在 残 差 图 像 中 当前 — 据 可 以 看 出 : 
像素 点 坐标 为 j) ,将 其 与 间隔 (m,n) 位 置 的 像素 点 亮度 值 求 DEEH DM ,的 变化 非常 有 规律 。 在 迭代 过 程 中 的 初始 阶 
积 ， 重 复 对 ;中 所 有 像素 点 做 上 述 运算 并 求 和 ， 可 以 表示 为 段 ， 度 量 值 变化 比较 大 且 下 降 非 常 快 ， 这 说 明 图 像 质量 提升 很 
快 ， 随 着 迭代 次 数 的 增加 ， 度 量 值 pm 下 降 得 越 来 越 慢 ， 说 明 


Coe, (m,n) = 22450 Dr-m,i-n) 


(13) 


图 


像 质量 趋 于 稳定 。 
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b) 对 于 不 同类 型 模糊 、 不 同 模糊 程度 和 不 同 图 像 内 容 的 


图 


像 ， 很 难 依据 值 DM,, 本 身 大 小 来 决定 是 否 迭 代 过 程 应 该 终止 。 
例如 ， 表 2 中 度量 值 DM ,在 0.38、0.05 和 0.59 处 趋 于 稳定 ， 


ive EBRTHI 
aA VEL (FER, 


: Chin j 
江 顺 亮 ， 等 : 适用 于 和 迭代 型 去 模糊 算法 的 自 适应 先 


因此 无 法 直接 根据 度量 值 DM,, 的 大 小 决定 迭代 过 程 是 否 应 该 


终止。 


表 2 均匀 模糊 (K=[5 5]) 下 部 分 迭代 步 数 的 度量 值 DM 变化 情况 


Step 1 5 10 20 50 120 240 360 480 720 
DM, 5.588 0.92 023 0.29 035 0.36 038 038 038 038 
43 高 斯 模糊 (K=[15 15], S-3) 下 部 分 迭代 步 数 的 度量 值 PDM, 变化 情况 
Step 1 5 10 20 50 120 240 360 480 720 
DM, 2703 0.02 0.03 0.04 0.04 0.04 0.05 0.05 0.05 0.05 
表 4 运动 模糊 (A=0, L=10) 下 部 分 迭代 步 数 的 度量 值 DM ,. 变化 情况 
Step 1 5 10 20 50 120 240 360 480 720 
DM, 1410 0.78 0.59 0.49 0.48 0.54 0.59 0.59 0.59 0.59 


2.3 自 适应 迭代 终止 条 件 


从 上 文 给 出 的 去 模糊 度量 DM ,的 值 随和 迭代 次 数 变化 的 趋 


势 来 看 ， 很 难 依据 评价 度量 值 本 身 值 的 大 小 来 决定 是 否 迭 代 过 


程 


应 该 中 止 。 为 此 ， 


本 文 提出 依据 DM,, 值 变化 趋势 确定 AISC 


条 件 ， 即 通过 检测 连续 步 数 Step 范围 内 DM ,, 值 变化 的 绝对 值 


是 否 


否 小 于 预 设 的 闵 值 。 来 判定 
AISC 条 件 满足 ， 则 去 模糊 算法 的 迭代 过 程 终止 。 基 
数 Step 和 净值 e 是 AISC 条 件 两 个 重要 的 参数 ， 它 们 


图 像 质量 是 否 已 经 趋 于 稳定 。 若 
此 ， 连 续 步 


体 值 的 


设置 主要 根据 特定 去 模糊 算法 的 欠 代 收敛 速度 ， 依 据 大 量 实验 


为 了 设置 一 个 合理 的 步 数 和 闵 值 ， 选 取 如 图 3 中 10 张 图 


像 进行 实验 分 析 ， 图像 均 来 自 目前 文献 中 常用 的 TID2013 数据 


图 3 TID2013 图 


LIMES GUI 


有 模糊 (K=[9 9],S=3) 在 连续 搜索 步 数 Step 为 


15. 20. 25 下 进行 ， 浆 值 e 的 值 设 置 为 0.01、0.02 和 0.03。 表 


5 数据 为 10 张 图 像 终止 时 与 720 次 终止 时 
给 出 了 PSNR, SSIM, FSIM 


均值 ， 表 5 


图 像 质量 值 差 值 的 平 


值 三 种 图 像 质 量 评价 


指标 的 差 值 ， 其 中 正 值 表示 采用 本 文 所 提出 的 AISC 条 件 后 所 
获得 的 图 像 质量 优 于 NCSR 算法 原 720 步 固 定 迭 代 次 数 所 获得 


的 结果 ， 负 


WE e 被 设置 为 0.01 时 ，FSIM 值 才 
为 20 时 在 e =0.01 时 最 好 的 去 模糊 效 


为 正 


值 则 表示 要 差 。 表 5 中 加 粗 加 下 画 线 的 为 PSNR、 
SSIM, FSM 值 在 不 同 设置 条 件 下 的 最 佳 


E. ERKE, RA 
值 ， 且 连续 步 数 Step 


H 
ES 


忆 此 本 文选 用 步 数 为 


20. REX 0.01 作为 迭代 终止 条 件 的 参数 。 下 文 所 有 的 实验 数 


据 都 是 在 这 个 参数 设置 条 件 下 获得 的 


像 集 [20] 中 10 张 图 像 


o 


de 5 AERP T EL NDS TNR LESE SESS Se BL 28. 720 步 质 量 值 的 平均 误差 
BU [Ei 0.01 0.02 0.03 
步 数 PSNR SSIM FSIM PSNR SSIM FSIM PSNR SSIM FSIM 
15 02817 0.0407 0.0001 0.2456 0.0500 -0.0012 0.3111 0.0550 -0.0008 
20 0.2964 0.0373 0.0017 0.1969 0.0469 -0.0020 0.2501 0.0516  -0.0019 
25 0.2670 0.0343 0.0010 0.1659 0.0449 -0.0025 0.2139 0.0493  -0.0025 
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3 ”实验 结果 及 其 分 析 


3.1. 实验 测试 平台 和 运行 环境 
为 了 验证 方法 的 有 效 性 和 稳定 性 ， 选 取 两 个 具有 代表 性 的 


Couple, Elk, Goldhill, House, Lena, Man 和 Peppers 图 像 ; 
第 二 个 图 像 集 由 上 文 使 用 到 的 TID2013P9 中 的 10 张 图 像 和 文 


3.2 ”常用 图 像 集 


为 了 验证 AISC 和 迭代 终止 条 件 的 有 效 性 ， 首 先 在 常 / 
集 上 进行 测试 ， 测 试 结果 如 表 6-8 所 示 。 从 表 6~8 中 数据 可 
看 出 : a) 终 止步 数 分 布 在 200 步 以 内 ， 远 少 于 原 NCSR 算法 
720 次 固定 的 迭代 次 数 ; b) 终 止步 数 对 应 的 


步 质 量 值 ( 即 原 NCSR 算法 获得 的 图 像 质量 值 ) 差 别 不 大 ， 人 


肉眼 其 实 比 较 难于 区 分 。 


从 左 到 右 分 别 对 应 的 是 Elk 


体 的 分 为 两 种 情况 ， 如 图 6 所 示 


图 像 在 参数 为 (K=[15 15],S=3) 的 


图 像 质量 值 与 第 720 


图 像 集 上 进行 测试 : 第 一 个 图 像 集 由 各 大 文献 中 广泛 使 用 的 测 
试图 像 构成 ,如 图 4 所 示 , 这 些 图 像 是 Barbara, Boat, Cameraman, 
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献 [17] 中 的 90 张 图 像 (部 分 图 像 如 图 5 所 示 ) 构 成 。 由 于 图 像 集 
中 内 容 复 杂 、 纹 理 特征 丰富 ， 本 文 称 之 为 挑战 性 图 像 集合 ， 非 
常 适合 用 来 验证 所 提出 的 AISC 条 件 的 鲁 棒 性 。 算 法 运行 的 硬 
件 平 台 为 Intel Core™ i7-6700 3.40 GHz CPU 处 理 器 ，8 GMB 
内 存 。 软 件 系统 配置 为 Windows? 操作 系统 ， 编 程 环 境 为 
MATLAB 2014a. 为 了 说 明 迭 代 终 止 条 件 的 有 效 性 ,采用 PSNR, 
SSIMISI 和 FSIM09 图 像 质量 评价 指标 衡量 去 模糊 后 图 像 质量 


图 5 挑战 性 图 像 集中 10 KARKI 


TENZ 


以 
的 


的 


, 


E 


In] 


"ERS ES 


斯 模糊 下 原始 未 失真 图 像 、 模 糊 图 像 、 自 适应 迭代 终止 步 数 对 
应 的 图 像 和 第 720 步 对 应 的 图 像 。720 PIKIR TER M ERA 
量 较 自 适应 终止 获得 的 图 像 质量 稍 有 改善 ， 但 不 易 被 觉察 到 ; 

如 图 7 所 示 ，Lena 图 像 参 数 (K=[5 5]) 的 均匀 模糊 下 原始 未 失真 
图 像 、 模 糊 图 像 、 自 适应 迭代 终止 步 数 对 应 的 图 像 和 第 720 步 
对 应 的 图 像 。720 步 迭 代 完 成 后 的 图 像 质 量 较 自 适应 终止 获得 
的 图 像 质量 稍 差 , 但 也 不 易 被 觉察 到 。 由 此 可 见 , 采用 AISC 终 
止 条 件 后 ， 改 进 算法 可 以 适时 终止 迭代 过 程 ， 极 大 提高 了 计算 
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效率 。 在 图 像 质量 方面 存在 提高 和 降低 两 种 可 能 ， 但 是 这 两 种 
情况 与 原 算法 所 获得 的 图 像 在 图 像 质量 方面 的 差别 很 小 ， 人 的 
肉眼 较 难 区 分 。 改 进 算法 在 不 损失 太 多 图 像 质量 的 前 提 下 实现 
了 大 幅度 提高 计算 效率 的 目的 。 
表 6 均匀 模糊 下 自 适应 迭代 终止 对 应 的 质量 值 与 

第 720 步 质量 值 的 平均 误差 比较 


模糊 类 型 中 止步 数 范 围 PSNR SSIM — FSIM 

均匀 模糊 (K=[5 5]) 41~54 0.25 0.0135 0.0030 
和 匀 模 糊 (K=[7 7]) 16~53 -0.42 -0.0032 -0.0085 (EE E 

均匀 模糊 (K=[9 9]) 12-47 -0.85  -0.0142  -0.0166 


7 ”高 斯 模糊 下 自 适应 迭代 终止 对 应 的 质量 值 与 


第 720 步 质量 值 的 平均 误差 比较 

模糊 类 型 中 止步 数 范围 PSNR SSIM FSIM H 

高 斯 模糊 (K=[5 5],S=3) 70~85 1.05 0.0427 0.0173 i 

高 斯 模糊 (K=[9 9],S=3) 6-16 -0.85 0.0401  -0.0248 j H 

高 斯 模糊 (K=[15 15],S=3) 6~15 -1.55 -0.0537 -0.0494 i 
(@) 和 迭代 终止 步 数 (第 41 步 ) 对 应 的 图 像 (OF 720 步 对 应 图 像 
和 图 7 均匀 模糊 (K=[5 5]) 条 件 下 Lena 图 像 的 视觉 效果 对 比 

第 720 步 质量 值 的 平均 误差 比较 
模糊 类 型 中 止步 数 范围 “PSNR SSIM — FSIM 3.3 ”挑战 性 图 像 集 

PR 为 了 进一步 验证 AISC 是 否 具 有 较 强 的 鲁 棒 性 ， 在 挑战 性 


ZARHAL 。 91150 — 027 00] ooo 图 像 集合 上 完成 了 实验 。 表 9-11 列 出 的 是 三 种 模糊 类 型 下 的 
IEXUMMHA-OL-4) — 77-90 005 00090 00020  ” 实验 结果 ， 整 体 情况 与 常用 图 像 集 的 情况 类 似 。 总 体 来 讲 ， 本 
文 提出 的 AISC 迭代 终止 条 件 能 够 在 保证 图 像 质量 损失 较 少 的 
情况 下 (也 有 部 分 情况 提升 图 像 质 量 )， 对 算法 迭代 适时 的 终止 ， 
有 效 地 改善 了 计算 效率 。 


表 9 均匀 模糊 下 自 适应 欠 代 终止 对 应 的 质量 值 与 
第 720 步 质量 值 的 平均 误差 比较 


模糊 类 型 中 止步 数 范围 PSNR SSIM FSIM 
均匀 模糊 (K=[5 5]) 39~63 0.22 0.0161 0.0208 
均匀 模糊 (K=[7 7]) 37~60 -0.32 0.0012 -0.0058 
均匀 模糊 (K=[9 9] 15-55 -0.72 -0.0070 -0.0133 


(b) 模 糊 图 像 


表 10 高 斯 模糊 下 自 适 应 迭代 终止 对 应 的 质量 值 与 
第 720 步 质量 值 的 平均 误差 比较 


模糊 类 型 中 止步 数 范围 PSNR SSIM — FSIM 
高 斯 模糊 (K=[5 5],S=3) 57~91 1.01 0.0441 0.0175 
高 斯 模糊 (K=[9 9],S=3) 4-36 -0.47 0.0556 -0.0080 
高 斯 模糊 (K=[15 15],S=3) 4~19 -1.33 -0.0458 -0.0442 
(c) 和 迭代 终止 步 数 (第 10 步 ) (d) 第 720 步 对 应 图 像 
对 应 的 图 像 


图 6 高 斯 模糊 (K=[15 15],S=3) 条 件 下 Elk 图 像 的 视觉 效果 对 比 


录用 稿 
表 11 运动 模糊 下 自 适应 迭代 终止 对 应 的 质量 值 与 
第 720 步 质量 值 的 平均 误差 比较 
模糊 类 型 中 止步 数 范围 PSNR SSIM FSM 
运动 模糊 (A=0,L=10) 67~164 0.58 0.0234 0.0091 


运动 模糊 (A=0,L=20) 22~1 


51 


运动 模糊 (A=0,L=40) 22~76 


0.25 | 0.0206 0.0076 
-0.64 0.0213 -0.0075 


3.4 执行 时 间 


为 了 测试 引入 去 模糊 效果 度量 后 所 增加 的 运行 时 间 ， 以 一 


其 最 大 的 计算 量 在 于 中 间 估 计 


故 执行 时 间 不 长 。 而 现 有 的 评 


张大 小 256*256 的 Lena Ef ill 
执行 一 步 和 计算 去 模糊 度量 所 需 执行 时 间 , 列 在 表 12 中 。 由 表 
12 中 可 知 , 本 文 所 提出 的 去 模糊 图 像 质 量度 量 所 增加 的 时 i 
常 少 (都 在 毫秒 量 级 上 ), 几乎 可 以 忽略 不 计 。 其 原因 在 于 本 文 所 
提出 的 去 模糊 效果 度量 是 专门 针对 去 模糊 算法 的 特点 而 设计 的 
图 像 与 模糊 核 的 卷 积 操作 ， 而 卷 
积 操 作对 于 目前 的 计算 机 来 说 有 大 和 
介 准确 性 比较 好 的 NR-IQA 算法 


试图 像 ， 记 录 NCSR 内 循环 


| 


非 


量 成 熟 的 加 速算 法 可 以 利 


很 多 运行 时 间 都 在 秒 的 量 级 上 (有 的 


其 至 高 达 几 十 秒 ), 因此 , 本 


文 所 提出 的 DM 度量 远 比 引入 


模糊 图 像 质量 还 有 一 定 程 度 上 


升 。 


BUG If] NR-IQA 算法 所 增加 的 执 
行 时 间 要 少 很 多 。 另 外 由 上 文 可 知 ， 
算法 一 般 情况 下 在 迭代 200 步 以 内 就 会 终止 ,因此 执行 时 间 大 
为 缩短 。 其 实 ， 不 仅仅 是 在 内 循环 中 减少 了 和 迭代 步 数 所 获得 的 
提升 , 在 外 循环 的 中 所 节省 的 字典 


采用 ASIC 条 件 后 NCSR 


学 习 (dictionary learning) 时 间 


也 是 非常 可 观 的 。 因 此 , 引入 AISC 条 件 后 , NCSR 算法 在 计算 
性 能 方面 的 提升 是 非常 可 观 的 ， 而 且 在 多 数 时候 ， 所 获得 的 去 
虽然 在 茶 些 时 候 所 获得 的 图 


像 质 量 不 如 原 NCSR 算法 ， 但 在 视觉 上 差别 也 不 大 。 


表 12 NCSR 算法 和 从 代 一 步 所 需要 时 间 和 
计算 DM 度量 所 需 时 间 对 比 /ms 


Blur type Parameters 


One step DM metric 


Uniform blur K=[5 5] 


Gaussian blur K-[9 9], S=3 


Motion blur A=0, L=40 


102.1 4.4 
88.7 4.3 
69.3 4.4 


4 ”结束 语 


现 有 基于 迭代 求解 方式 的 去 模糊 算法 在 实现 时 采用 固定 的 


迁 代 次 数 ， 并 不 能 保证 所 获得 


存在 着 算法 的 时 间 效 率 比 较 低 的 问题 。 本 文通 过 分 析 残 差 
的 特点 ， 利 用 残 差 图 像 中 各 个 像素 点 之 间 的 亮度 值 差 异 作为 衡 


的 去 模糊 图 像 质量 是 最 佳 的 ， 且 


ERI 


像 


量 去 模糊 效果 的 度量 值 ， 并 基 


lm 


自 适应 地 终止 迭代 过 程 ， 改 善 


于 这 个 度量 设计 了 ASIC 自 适 应 
的 迭代 终止 条 件 。 采 用 ASIC 迭代 中 止 条 件 后 的 NCSR 算法 较 


原来 采用 固定 迭代 次 数 的 方法 在 保证 图 像 质量 的 前 提 下 ， 能 够 


$ 


了 原 算法 耗 时 长 的 问题 。 需 要 特 
别 说 明 的 是 ， 本 文 所 提出 的 ASIC 和 迭代 中 止 条 件 并 非 只 能 用 于 


NCSR 算法 ， 其 工作 原理 适 | 


j 于 所 有 和 迭代 型 的 去 模糊 算法 。 当 


N 


, 


i* 
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然 它 的 两 个 关键 参数 需要 根据 茶 个 特定 去 模糊 算法 在 迭代 过 程 
中 的 收敛 速度 合理 设置 后 方 可 应 用 。 
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